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Abstract 

A number of extensions to the Standard Model of particle physics predict a permanent electric dipole moment of the electron 
(eEDM) in the range of the current experimental limits. Trapped ThE^ will be used in a forthcoming generation of the JILA eEDM 
experiment. Here, we present extensive survey spectroscopy of ThE^ in the 700 - 1000 nm spectral region, with the 700 - 900 nm 
range fully covered using frequency comb velocity modulation spectroscopy. We have determined that the ThE^ electronic ground 
state is X ^Ai, which is the eEDM-sensitive state. In addition, we report high-precision rotational and vibrational constants for 14 
ThE^ electronic states, including excited states that can be used to transfer and readout population in the eEDM experiment. 
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1. Introduction 

A measurement of the permanent electric dipole moment of 
the electron (eEDM) is a direct test of parity- and time-reversal 
symmetry violation iniEiii. Electron EDM experiments have 
typically been performed on beams of neutral atoms il or neu¬ 
tral molecules 13 a, with the current experimental limit of 
\de\ < 9.7 X 10“^® e-cm set by the ACME collaboration’s molec¬ 
ular beam experiment IT]. The JILA eEDM experiment in¬ 
stead uses trapped molecular ions, leveraging the large effective 
electric held IS] [51 and suppression of eEDM systematic errors 
available in molecular ^Ai states nniiiii, while also achieving 
a long measurement coherence time IfT^ . The JILA experiment 
is currently using trapped HfF^ ca, but ThE^ is a promising 
candidate molecule for a future generation experiment due to 
its larger effective electric held Ha (13 da and longer-lived 
^Ai state. Assigning the ground state of ThE^ has been exper¬ 
imentally challenging due to the proximity of both and 
states inidii- Calculations suggested that ^Ai might be the 
ground state but uncertainties prevented a conclusive determi¬ 
nation 113 dS- Here, we show that the ThF^ ground state is 
X ^Ai, and that the a state lies 314 cm“' higher in energy. 
This is benehcial for the eEDM experiment, since the measure¬ 
ment coherence time will not be limited by the spontaneous 
lifetime of the “science” state. In addition to ^Ai and *2^, we 
have measured high-precision rotational and vibrational con¬ 
stants for 12 additional electronic states. The results reported 
here are benchmarks for ab initio theory, and may assist in un¬ 
derstanding actinide bonding. 
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The lack of spectroscopic information about electronic states 
of ThF'^ above 10,000 cm ', and the difficulty of perform¬ 
ing survey spectroscopy on a small number of trapped ions, 
necessitated a tool that provides broad bandwidth, high pre¬ 
cision, high sensitivity, and ion-selectivity. Frequency-comb 
velocity-modulation spectroscopy (comb-vms) ll20]l2Tll . which 
combines cavity-enhanced direct frequency comb spectroscopy 
1231231 and velocity-modulation spectroscopy (vms) llllES, 
fulhlls these requirements. With this system we have scanned 
the entire wavelength range of 696 - 900 nm, acquiring 3,260 
cm“* of continuous spectra at 150 MHz (0.005 cm“') step-size. 
Additionally, guided by the results of the comb-vms and the 
PFI-ZEKE spectra from Barker et al. IfTTl . we have performed 
targeted cw-laser vms scans to identify additional experimen¬ 
tally useful transitions and to assign electronic-state Q values 
via observation of low-J rotational lines. In total, we have fit 
27 rovibronic transitions in ThF'^ and 5 transitions in a differ¬ 
ent thoriated molecular ion, potentially ThO^. We discuss the 
production and detection of ThF^, details of the ground state 
assignment, and the assignment of other ThF^ electronic states. 

2. Experimental details 

2.1. Frequency comb vms 

Fig. 0 shows an overview of the comb-vms system. A 3- 
GHz repetition-rate titanium-doped sapphire frequency comb 
is spectrally broadened in nonlinear hber and coupled into an 
enhancement cavity surrounding a 1-kW electrical discharge 
cell. A VIPA imaging system resolves individual 

comb modes and provides 1,500 simultaneously resolved comb 
modes per image. 

The spectroscopically usable fundamental bandwidth of the 
comb is 780-850 nm, centered at 815 nm. To obtain maximum 


Preprint submitted to 


September 15, 2015 






f 

rep 


Figure 1: Experimental schematic. A 3 GHz Ti:Sapph comb is broadened 
in nonlinear fiber and coupled into a bow-tie ring cavity, propagating in either 
direction through a discharge of ThF 4 . Four liquid crystal variable retarders 
and a polarizing beam splitter control the direction of light propagation through 
the discharge. The cavity free spectral range is 1/25 of the comb repetition 
frequency. After exiting the cavity, light is sent to a cross-dispersive VIPA 
imaging system with single comb mode resolution and imaged onto a Heliotis 
lock-in camera da Approximately 1,500 comb modes are resolved in a single 
image. 


power throughput to the experiment, the output of the frequency 
comb is immediately coupled into the nonlinear fiber to broaden 
the spectrum. A half-wave plate before the nonlinear fiber in¬ 
put enables control of the output spectrum ||28]| . No single fiber 
gave a smooth, high power spectrum in the entire 700 - 900 
nm range, so we used three different fibers, replacing them as 
needed to cover the entire region. A Newport SCG-800-CARS 
fiber resulted in significant power in the 850 - 900 nm spec¬ 
tral region but practically no light in the 700 - 780 nm region. 
85 cm of Thorlabs NL-2.8-850-02 fiber, pumped in the normal 
dispersion regime, provided a small amount of broadening and 
resulted in high power in the 730 - 780 nm range. 80 cm of 
Thorlabs NL-2.4-800 fiber, pumped close to the zero dispersion 
wavelength, resulted in most of the laser power being shifted to 
the 680 - 730 nm regime. 

The fiber output spectrum often contained a significant frac¬ 
tion of laser power shifted to a spectral region outside of the 
desired scanning range. After the fiber, low- and high-pass op¬ 
tical filters were used to select a 50 nm spectral window. The 
spectrum after the filters was optimized and monitored using 
a grating spectrometer with a USB readout. A linear polarizer 
was placed after the filters to eliminate the different polariza¬ 
tion components due to nonlinear polarization rotation in the 


fiber. An electro-optic modulator, driven at 17 MHz, applies 
RF sidebands to each comb tooth. 

After the modulator, the comb light is fiber-coupled to the 
table supporting the discharge cell. A 2.5-meter-long low- 
finesse (F ^ 70-100) enhancement cavity surrounds the dis¬ 
charge tube. Every 25th cavity mode is coupled to a comb 
mode, and the cavity is locked to the comb using the Pound- 
Drever-Hall technique Il29l . In order to obtain maximum cav¬ 
ity transmission, AR-coated windows were placed on each end 
of the discharge apparatus (R < 0.1% per surface 700 - 900 
nm), and T = 3% cavity input couplers were used, running 
the cavity in the over-coupled regime. Liquid-crystal variable 
retarders and a polarizing beam splitter were used to rapidly 
switch the direction of light propagation through the cavity and 
the discharge tube. By subtracting the spectrum acquired from 
each direction, common-mode laser noise and asymmetrical 
discharge effects, such as an asymmetrical concentration mod¬ 
ulation of molecules, are suppressed. 

After exiting the cavity, the comb light is imaged using a 2D 
VIPA spectrometer Eol with single comb mode resolution onto 
a (Heliotis) lock-in camera ll26l . To fully resolve different VIPA 
orders, a 1800 lines/mm diffraction grating was used for the 780 
- 900 nm spectral range, and a 2160 lines/mm grating was used 
for the 700 - 780 nm range. A cw titanium-doped sapphire laser 
is co-propagated with the comb light for frequency calibra¬ 
tion. The cw laser is locked to its internal reference cavity, and 
the frequency is measured by a 100-MHz-accuracy wavemeter. 
With interleaving scans, we obtain a maximum of 150 cm“' of 
continuous spectra over 1,500 simultaneous spectral channels at 
150 MHz (0.005 cm“') step-size in a 25 minute total data col¬ 
lection time, with 3 x 10“^ Hz“''^^ (spectral element)^''^^ frac¬ 
tional absorption sensitivity. The acquisition time is currently 
limited by software and with improvements could be shortened 
to 10 minutes. 

The discharge tube is a straight 50 cm length of 1.7 cm inner- 
diameter, 1.9 cm outer-diameter quartz tubing. The tube is cen¬ 
tered inside a home-built oven, which is heated to 920 °C to 
obtain a sufficient vapor pressure of ThF 4 ll^ [32l . Initially, 
smaller diameters of quartz tubing were used but were rapidly 
blocked during operation. The blockage is likely due to the re¬ 
action of ThF 4 with quartz at high temperature to form ThOa 
1331, which has a low vapor pressure at 920 °C lf34l and should 
rapidly plate out along the discharge tube. A straight alumina 
(AI 2 O 3 ) tube was tested but the discharge ceased to run down 
the length of the tube above 800 °C, instead opting to discharge 
to the grounded gas lines. 

For each oven run, approximately Ig of ThF 4 is loaded into 
the discharge tube, and the tube is pumped out with a scroll 
pump to 10 mTorr. Helium gas flows through the tube at a pres¬ 
sure of 5 Torr. Three heaters are used to maintain a temperature 
of 920 °C in the center of the oven and a temperature of 880 
°C on either edge of the oven. Higher temperatures were found 
to deplete the ThF 4 faster but not to yield a significant increase 
in ThF^ signal. Single oven runtimes were typically 3-4 hours 
before the signal fell to 80% of its peak value due to depletion 
ofThF4. 

During some scans we observed bands that we tentatively 


2 




















attribute to ThO^ (see Section 3.31. Attempts were made to 
remove the ThO^ signal from the discharge while monitoring 
ThF^ and ThO^ absorption strengths with cw-vms (see next 
section). Running the discharge with different buffer gases, 
such as Ne or Ar, resulted in the loss of all ThF^ signal while 
retaining the same ThO^ signal. Contrary to our expectations, 
seeding the discharge with SFg gas again resulted in the loss 
of all ThF^ signal and doubled the ThO^ line strengths. In all 
cases switching back to a He discharge resulted in an immedi¬ 
ate return to prior ThF^ line strengths. The ThF^ signal was 
only positively affected by total oven runtime, 

i.e. time when the oven was at 920 °C with ThF 4 loaded 
and a discharge running. After roughly 5 total hours of oven 
runtime, the ThO^ signal decreased by a factor of two and the 
ThF^ signal increased by a factor of two. This is likely due 
to Th-based discharge byproducts coating the tube walls and 
suppressing the ThF 4 - quartz reaction. 


2.2. cw Ti. Sapph and diode laser vms 

Twenty-one vibrational bands were ht in the comb-vms spec¬ 
trum. The total spectral coverage of the comb-vms system is 
limited by coatings on the VIPA etalon and by the decreased 
quantum efhciency of the silicon-based lock-in camera at wave¬ 
lengths longer than 900 nm. Cw-laser vms was used to extend 
spectral coverage beyond the 700 - 900 nm range in order to 
identify additional favorable electronic states for coherent pop¬ 
ulation transfer between X-a and to increase the signal-to-noise 
ratio on low-J rotational lines. For example, low-J lines of the 
Q = 0“ <—^A] (0, 0) band were observed at 685.4 nm, a wave¬ 
length inaccessible with the current optics coatings in the comb- 
vms system (see Section [XT] ). 

Due to the lower complexity of cw scans, a higher single¬ 
pass signal-to-noise ratio can be obtained. However, the total 
data acquisition time is more than 20 times slower per band¬ 
width than comb-vms. For the cw scans, the enhancement cav¬ 
ity was removed, and the laser was split into two components 
with a power ratio of approximately 2:1 that counterpropagated 
through the discharge tube. Roughly 1 mW of power is con¬ 
tained in the lower power beam. The lock-in signal is detected 
on the difference port of an autobalancing detector llTSl . The 
best scan rate achieved is 1 cm“' in 4 minutes with a 300 ms 
lock-in time constant, limited by reading the wavemeter and 
scanning slowly enough to maintain laser lock. The best single¬ 
pass fractional sensitivity reached is 2 x 10 “^ about a 

factor of 2 above shot noise at 1 mW, limited by a combination 
of electrical and acoustic pickup from the discharge. 

Low-J lines of the Q = 0 <— ^Ai (0, 0) band at 685.4 nm 
were scanned using an external-cavity diode-laser with a 680 
nm diode. The high-precision rotational and vibrational con¬ 
stants measured from the Q = 0“ <— ^Ai (0, 1), (1, 2), (2, 3), 
(0, 2), and (1, 3) bands allowed us to pinpoint line centers in 
the (0, 0) vibrational band to 0.1 cm ' accuracy. Feed forward 
to the diode current resulted in a mode-hop free scan range of 
0.2 cm"'. Fractional sensitivity with 50 pW optical power was 
limited by electrical pickup to 1 x 10 "^ Hz"*^^, yielding signal- 
to-noise > 10 on the weakest low-J lines with a 1 s lock-in time 
constant. 


3. Results and analysis 

We have observed 14 different ThF^ electronic states. To in¬ 
dicate state assignments Hund’s case (a) notation is used when 
possible, and case (c) notation is used when only the Q-value is 
known. All of the acquired data is shown in Fig. |^(a); a total 
of 3,260 cm"' was acquired with comb-vms (black), and 400 
cm"' with cw-vms (red). ThF^ bands often presented with at 
least a single P-, Q-, and R-branch, and frequently had two of 
each branch due to A-doubling in both electronic states. Fig. 
|^(b) shows a typical six branch vibrational band, with a clear 
R-branch band head and a strong Q-branch. The A-doubling is 
strong enough to fully resolve the two parity components by J 
= 3. 

For an observed vibrational band, rotational constants are de¬ 
termined by htting assigned line centers, dehned as the zero¬ 
crossing, to a general energy expression: 

E^vo + (B' - s'k'l2)J'(J' -Hi) - (D' - s'k'ol2)J'^{J' + \f 
- [{B" - s''k" I2)J"{J" -Hi) - {D" - s"k'^l2)J"Hj'' + 1)2]. 

( 1 ) 

Here, i is a parity term: i = -Hi for e-symmetry and i = -1 for 
/-symmetry, and k and k^ are generic A-doubling constants pro¬ 
portional to J(J-h 1) and J2 (J-h 1)2, respectively. Bands with A- 
doubling observed in each electronic state have all six branches 
ht simultaneously, with s = +1 corresponding to the two parity 
levels in each state. Most observed transitions have the promi¬ 
nent A-doubling contribution arising from the upper state; in 
these transitions k" and k" are set to 0 . 

Vibrational constants are determined from the relationship. 


v.-.. = r; - T'J + co'^iv' + 1/2)- aj,x'e(y' + 1 /2)' 

- [w"(v' + 1/2) - aj,Xe(y' + 1/2)"], 

where Vy/,,// is the (v', v") band origin, ate is the vibrational con¬ 
stant, ateXe is the anharmonicity constant, and Tg is the energy 
to the minimum of the electronic potential. For states with too 
few transitions observed to directly determine ate and ateXe, a 
Morse potential was assumed, allowing the substitution 


tOeXe 


Q'gtUg CTeOte 

3^e ^ 


+ Be, 


(3) 


where Og is the rotational-vibrational coupling constant and Bg 
is the equilibrium rotational constant. 

Vibrational assignments are determined by comparing ob¬ 
served band intensities to Franck-Condon factors and a thermal 
vibrational distribution. With prior knowledge of vibrational 
constants from Barker et al. HU, and our measured B val¬ 
ues, Franck-Condon factors were calculated and compared to 
observed intensities for the initial bands measured. With the 
rotational and vibrational constants we have observed in ThF^, 
the vibrational (0, 0) is always the strongest band. The highest 
intensity bands in the Av = +1 vibrational series are about 1.5-2 
times weaker than the ( 0 , 0 ) band. 

Due to the high density of electronic states, the relatively 
small difference in vibrational constants for different electronic 
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Figure 2: All acquired ThF^ spectrum. The energy level diagram on the right indicates all experimentally observed electronic energy levels. All y-axes are in 
10“^ single pass fractional absorption. All x-axes are in cm“*. (a) All ThF^ spectra acquired with comb-vms (black) and cw-vms (red), (b) Fitted positions to 
III [12.99] «— ^A 2 (0, 0). Note the presence of six branches — two P- (purple and green), two Q- (dark blue and light blue) and two R- (red and orange) - due to 
A-doubling in each electronic state, (c) Low-J lines in the = 0“ ^ Ai (0, 0) band. The solid colors indicate different branches for an £2' = 0 «— £2" = 1 transition 

(P- green, Q- blue, R- red). The dashed orange lines represent the R-branch intensity if the transition was instead £2' = 1 <— £2" = 0. In each case the overall band 
intensity is set by scaling the line intensities to match the strongest R-branch lines, which are off-scale peaks marked by solid red lines. P(l) is present, and R(0) is 
clearly absent. 


states, and the high oven temperature, the ThF^ spectrum is of¬ 
ten very congested. Fits are performed by manually identifying 
branches and assigning line centers to the rotational structure 
given by Eqn. 

We have developed a simple automatic line finder to assist 
with identifying bands in the most congested regions. Local 
maxima are identified and then matched with the nearest zero 
crossing of the correct line-shape phase within the observed 
linewidth of 1 GHz. A minimum intensity threshold is specified 
to avoid assigning noise to real lines, and a maximum intensity 
threshold is used to avoid assigning strong atomic features to 
molecular rotational lines. 

The identified lines are plotted in a Loomis-Wood style fig¬ 
ure, which displays J = 0-100 residuals on a single figure. Over 
such a large scale, pattern-forming rotational sequences can be 
easier to identify. By modifying B values, D values, and the 
band origin, vq, it is possible to precisely identify the line cen¬ 
ters of hundreds of rotational lines at once. Additionally, the 
large J-range of the plots enables clear identification of smaller 
high-J effects, such as A-doubling scaling as T^, which typically 
manifests in the spectrum as lines splitting into doublets above 


J = 40. 

3.1. ^Ai ground state assignment 

Assignment of the low-lying state is of paramount impor¬ 
tance to the eLDM experiment. Initial PLI-ZLKL experiments 
from Barker et al. El observed both the and ^ Ai states and 
measured the electronic energy difference to be 315(1) cm“^ 
However, due to the inability to fully resolve low-J lines, com¬ 
pounded by the similarity of B and aig values for the X and a 
states, it was not possible to conclusively assign either as the 
ThL^ ground state IfTTl . Our assignment necessitates observa¬ 
tion of an excited electronic state with coupling to both the 
and ^Ai states. Such a state will have a spectral signature of 
two vibrational bands with band origins separated by the X-a 
energy difference. In order to couple to both lower states, the 
upper state must be either Q = 0^ or Q = 1. The final assign¬ 
ment uses measurement of the A-doubling and an observation 
of low-J lines. 

In the comb-vms data we have fit transitions at vq = 
13935.996(7) cm“' and vq = 13286.716(7) cm“' with one set 
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Figure 3: Spectrum of ^Ai and coupling to a common H = 0^ state, (a) An = 0^ *— ^Ai vibrational band is observed at vq = 10471.889(7) cm“' (red 
axes). A strong Q-branch (blue markers) is observed in addition to the P- (green markers) and R-branches (red markers), (b) A transition is observed 314 cm“* 
lower in energy at vq = 10157.607(7) cm“* (blue axes). This transition has no Q-branch, indicating £1" = Q' = 0. The transition ordering indicates that ^Ai is 
lower in energy as illustrated in the energy level diagram on the right. 


of P-, Q-, and R-branches and a nonzero k constant. The pres¬ 
ence of a nonzero k and one of each branch implies that the two 
electronic states in the transition must be Q = 0 and Q = 1. The 
small magnitude of k, |A:| = 0.869(7)xl0“"^ cm“', makes it im¬ 
possible to differentiate between k in the upper state, k', or in the 
lower state, k". Because tOg for the and ^Ai states is so simi¬ 
lar it could not be used for electronic state identification. Based 
on calculated Franck-Condon factors and the intensity of the 
three transitions observed at vq = 13935.996(7), 13876.539(8), 
and 13817.279(8) cm“', and the two transitions observed at vq 
= 13286.716(7) and 13231.099(8) cm“', the transition at vq = 
13935.996(7) cm“' was identified as the Q = 0^ <— ^Ai (0, 1) 
vibrational band. 

The lower state rotational constant, Bq - 0.24264(3) cm“', 
and the vibrational constant, cOe - 656.96(1) cm both agree 
with Barker et al. im for low-lying states of ThF^. Also, the 

A-doubling proportionality constant, \k\ = 0.869(7) x 10“"^ 
cm *, is consistent with expectations for a An=i state: fi scaling 
and small magnitude 1^ . Thus, we tentatively assign the lower 
state to ^Ai. 

To confirm the lower state Q = 1 assignment we measured 
low-J lines of the Q = 0“ <—^Ai (0, 0) band at vq = 14589.09(2) 
cm ', as shown in Fig.[^(c). Because of the high oven temper¬ 
ature the spectrum is typically congested, but we were able to 
measure low-J lines of this band due to the high line strength. 
Higher-lying Q- and R-branch lines were scanned to confirm 
rotational constants and to fix the overall intensity of the band. 
Two compelling pieces of evidence are seen. First, the R(0) 
line, indicated by the dashed orange line, is absent, while the 
P(l) line, indicated by the solid green marker, is present. We 
clearly see R(l)-R(4). P(2)-P(4) are covered by strong Q- 

branch lines. Secondly, the Honl-London factors OtI for an 
Q" = 0 - = 1 transition considerably increase the intensity 

of the low-J R-branch lines relative to an Q" = 1 - Q' = 0 
transition, indicated by the dashed orange lines and the solid 
red lines, respectively. The observed R-branch line intensities 


agree with Q.” - 1. This confirms that the lower state is = 1 
and the upper state is Q = 0. 

We have observed another transition at vq = 10471.889(7) 
cm \ with one of each branch and nonzero k, shown in Fig. 
[^(a). The fitted values of B and D agree with Bq and Dq for 
^ 1 . Additionally, the magnitude of k agrees with the previ¬ 
ously measured transition at vq = 13935.996(7) cm“' but the 
sign is changed, which is expected if the upper state is of a dif¬ 
ferent overall parity, e.g. an Q = 0^ instead of an Q = 0“ state. 
Therefore, the lower state in this transition is also ^Ai. 

A vibrational band with no Q-branch is observed 314.282(7) 
cm * lower in energy, consistent with the X-a separation mea¬ 
sured by Barker et al. ifTTll . with vq — 10157.607(7) cm“'. The 
band origin is shown in Fig. [^(b). The absence of a Q-branch 
indicates the two electronic states in the transition both have Q. 
= 0, giving strong evidence for the lower state in this transition 
being To confirm the ground state assignment, the upper 
electronic state must be the same as the one observed in the vq 
= 10471.889(7) cm“* transition. 

With no Q-branch the band cannot be fit a priori. How¬ 
ever, by fixing B' and D' to the values measured in the vq = 
10471.889(7) cm“* transition, we obtained a good fit with rea¬ 
sonable residuals. The lower-state rotational constant fits to 
B” — 0.24542(3) cm“'. This is larger than the fit B(^Ai), as 
intuitively expected since an electron configuration (in 'S^) 
is typically more tightly bound than an sd configuration (also 
observed in HfF+ ||2T1 '). 

The fit residuals provide the final piece of compelling evi¬ 
dence for the upper electronic state being the same. Fit resid¬ 
uals for the Q = 0+[10.47] ^ (0, 0) at vq = 10471.889(7) 

cm“^ are plotted in Fig. |4](a). A strong avoided crossing per¬ 
turbation is evident at approximately J' - 62. Other ^Ai fits 
have no avoided crossing, so the perturbation must be present 
in the Q = 0^ upper electronic state. Fit residuals for the vq = 
10157.607(7) cm“' transition, shown in Fig. |^(b), display an 
avoided crossing of the same magnitude again at J' - 62. This 
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Figure 4: Fit residuals for ^Ai and coupling to a common upper state. 

The fit residuals of each transition exhibit an avoided crossing perturbation of 
the same magnitude at the same rotational quantum number, confirming the 
upper electronic state is the same in each transition. The dashed grey lines 
on each figure are guides to the eye and are identical, (a) Fit residuals for 
n = 0+[10.47]<-3 Ai (0, 0). (b) Fit residuals for n = 0+[10.47] <-(0,0). 

confirms that the upper state is the same for the two transitions, 
and thus that the assignments of and ^Ai are correct. Since 
the Q. - 0^[ 10.47] «— '2+ (0, 0) transition frequency is lower 
than that of the Q = 0^[10.47] <— (0, 0) transition, ^Ai is 

the electronic ground state. 

3.2. Additional ThF^ electronic states 

In addition to the and ^Ai electronic states we have ob¬ 
served rovibronic transitions that we have assigned to originate 
from the ^A2 and ^A3 states. These are given in Tableas the 
Hi [12.99] ^ ^A2, the Hi [14.22] ^ ^^2, and the ^U 2 <- ^As 
groups of transitions. 

For the ^A2 state, five vibrational bands of IIi [12.99] <— ^A2 
and one band of Hi [14.22] <— ^A2 were observed. Each band 
has two of each P-, Q-, and R-branch originating from a com¬ 
mon band origin, which indicates that A-doubling is present 
in each state. A small section of the Hi [12.99] <— ^A2 (0, 
0) band is shown in Fig. (b). Two strong Q-branches are 
seen, indicating that Q' + Q!'. The A-doubling has a strong 

dependence, consistent with the upper state being predomi¬ 
nantly nn=i Il38ll^ . B” - 0.24342(2), similar to the value of 
Be(^Ai) = 0.24311(7), and consistent with the measured value 
for ^A2 from Barker et al. lEl. Additionally, cj'J is consistent 


with the value measured by Barker et al. and again similar to 
W,(3Ai). 

To confirm this assignment, we observed low-J lines in the 
Hi [12.99] ^ ^A2 (0, 1), (0, 0), and (1, 0) bands as well as in 
the Hi [14.22] <— ^A2 (0, 0) band using the cw TitSapph. In each 
case, Q(l) is absent and Q(2) is present, thus leaving Q.” = 1 or 
2 as the only possibilities. R(l) is absent in each band, and the 
intensity of the low-J P- and R-branch lines is consistent with 
Q" > Q.'. Therefore, Q." = 2, Q' = 1, and the lower electronic 
state is ^A2. 

The assignment for is similar. Three vibrational bands 
of ^02 <— ^A3 were observed, with each band again having two 
of each P-, Q-, and R-branch at a common band origin. Here 
though the A-doubling has a strong dependence, as expected 
for a ^nn=2 upper state li3^ . B” - 0.24388(7) is similar to 
that of ^Ai and ^A2, and consistent with Barker et al. ini. 
Scanning near the band origins of the ^n2 <— ^A3 (0, 1) and 
(0, 0) bands with the cw laser revealed that P(2) is absent and 
P(3) is present. Low-J Q-branch lines could not be resolved due 
to contamination from other bands. However, the intensity of 
the low-J P- and R-branch lines is again consistent with fJ" > 
Q', concluding that Q." = 3 and Q' = 2. This confirms the 
lower state ^A3 assignment. The assignment of upper state ^n2 
is based on the strong, J'^ A-doubling, and the expected strong 
transition strength for a ^n2 <— ^A3 transition. 

Several bands exhibit strong perturbations at high-J, which 
manifest as fit residuals that deviate from Eqn. [^at 7 > 60. To 
obtain uniform fit residuals at high-J it was necessary to limit 
the range of J-values, or to include additional fit parameters in 
the Hamiltonian, such as HJ^(J + 1)^, or in some cases even 
higher powers of J. Higher vibrational levels of the Hi [12.99] 
state appear to become perturbed, likely as the result of prox¬ 
imity to a perturbing state near 14,000 cm“^. This is evident by 
three different observations. Eirst, k'j^ has a strong vibrational 
dependence, increasing from 2.1(2) x 10“^ cm“' in the v' = 1 
vibrational state, to 7.7(8) x 10“® cm“' in the v' — 2, and to 
10.1(7) X 10“^ cm“' in the v' - 3, the highest vibrational state 
we observe. Secondly, fit residuals deviate from Eqn. [T] at J 
> 80 for the v' = 2 and v' = 3 vibrational states, but not for 
the v' = 0 and v' = 1. Einally, assuming a Morse potential, 
cOe and cjeXe can be calculated for both the lower and upper 
electronic states using the (0, 1), (0, 0), and (1,0) bands. How¬ 
ever, the origins of the (2, 1) and (3, 2) bands are inconsistent 
with the values calculated, indicating that the upper state elec¬ 
tronic potential deviates from the Morse approximation beyond 
v' = 1. Similarly, is strongly vibrationally dependent in the 
^n2 <— ^A3 group of transitions, decreasing from 74(4) x 10“^ 
cm“* in the v' - 0, to 41(9) x 10“® cm“' in the v' - 1. This 
is likely due to a high density of unobserved (dark) electronic 
states. 

Perturbations are also evident in a number of other vibra¬ 
tional bands. The Q. = 0^ [14.29] <— ^Ai (0, 0) band shows a 
strong deviation from the expression of Eqn. above J - 50. 
The (0, 1) band has the same structure in the fit residuals. The 
(1,1) band, which ordinarily would be strong, is not seen, and 
the (2, 2) band is observed at vq = 13639.836(4) cm“', only a 
difference of 2 cm“^ from the band origin of the (0, 1). This 
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deviation from typical vibrational spacings indicates that the 
upper state electronic potential is heavily perturbed. Similarly, 
the III [14.22] <— ^A 2 (0, 0) band fit residuals deviate from the 
standard form above J - 50. Here, no other vibrational levels 
have been observed; the (0, 1) band is obscured by the strong 
Hi [12.99] <- (1, 0), and the (1, 1) and (2, 2) are not ob¬ 

served. 

Four additional ThF^ electronic states were observed but 
not conclusively assigned. In all cases, the lower electronic 
state B values are not consistent with the states observed by 
Barker et al. iflTl . This suggests that T^, > 4,000 cm“’, 
which in turn suggests that this population is strongly superther- 
mal. Franck-Condon factors indicate that the band at vq = 
10702.830(4) cm“* is a vibrational (0, 0) and the two bands 
at VO = 10085.683(7) and 10050.200(7) cm^' are the (0, 1) and 
(1, 2), respectively. A third band was observed at vq ~ 10,014 
cm *, consistent with the location and intensity of the (2, 3), 
but was not fit. Each transition has two of each P-, Q-, and 
R-branch, again indicating A-doubling in each electronic state. 
The lower electronic state exhibits /""^-scaling A-doubling with 
a strong vibrational dependence; k'^ - 34.0(7) x 10“® cm“' 
in the v” - 0, and it decreases to 14(1) x 10“® cm“^ in the 
v" - 1 and 5.5(3) x 10“® cm“' in the v" - 2. Low-J line in¬ 
tensities in the (0, 0) and (0, 1) bands reveal that Q" < Q'. 
R(l) is absent and R(2) is present. The magnitude and scal¬ 
ing of the A-doubling, the lower value of B, and the low- 
J lines suggest that the lower electronic state may be a ^n 2 
state. With the vibrational assignment we can calculate con¬ 
stants for the electronic states. For the lower state, we obtain 
B’’ = 0.23402(5) cm'^ w" = 619.8(2) cm-^, = 1-34(9) 

cm“*, and a" = 0.74(4) x 10“"^ cm“'. For the upper state 
we measure: R' = 0.22363(5) cm“*, w' = 593.1(4) cm“*, 
^eX'e - 1-70(9) cm“', and a' = 0.89(4) x 10“"^ cm“'. 

A pair of transitions was observed atvo = 12721.532(7) cm“' 
and Vo = 12661.053(7) cm '. Each band has a single P- and R- 
branch and a strong Q-branch. k and ko both fit to 0, indicating 
very small A-doubling in one or both of the states. Due to spec¬ 
tral congestion and the relatively low intensity of the band, no 
low-J lines were observed. The intensity of the bands is consis¬ 
tent with vibrational quantum numbers of Av = -H1 or Av = -1. 
A third band, likely the (v' + 2, v" -H 2), can be barely resolved 
at Vo ~ 12,600 cm“', but it is too weak to fit. The (0, 0) is 
typically the strongest band, but is not immediately observed 
when searching within a reasonable range of Wg values; it is 
likely obscured by the presence of stronger bands in these re¬ 
gions. Without a definitive vibrational assignment we cannot 
conclusively determine Rg, Wg, or cJeXe- 

3.3. ThO^ 

We have also fit five bands that do not originate from ThP^. 
Based on the rotational constants, these bands appear to be orig¬ 
inating from two different electronic states, separated by about 
582 cm \ coupling to a common upper state. The fit rotational 
constants are nearly 1.5 times larger than those for ThE^ values: 
R a! 0.34 cm ', compared to R(ThE^) ^ 0.24 cm"'. Addition¬ 
ally, the bands must be fit with half-integer J, indicative of a 
molecule with a spin doublet. ThO^, experimentally studied 


by Goncharov et al. in 2006 |[39l, seems a promising candi¬ 
date. However, our measured constants, given in Table do 
not form a fully consistent picture with the values measured by 
Goncharov et al., so we warn that our identification of ThO^ as 
the unknown species is tentative. 

ThO^ bands were fit to the form given in Eqn. except 
with J half-integer. All five transitions fit well out to J = 100 
with no A-doubling terms, i.e. k - ki, - 0. The vibrational 
assignment is tentative and is based on overall band intensities 
as predicted by Eranck-Condon factors. The electronic state 
with Ro = 0.33984(2) cm"' is consistent with the ^A 3/2 state of 
ThO^. R' = 0.32610(1) cm"' is a reasonable value for a higher- 
lying electronic state. However, the third electronic state, with 
Ro = 0.33900(2) cm"', is only consistent with a higher-lying 
or ^A 5/2 vibrational level. The upper electronic state in 
each group of transitions appears to be identical; R is consis¬ 
tent to within 0.00002 cm"'. If the upper state is the same, the 
two lower states must be separated in energy by 12146.812(1) 
cm"' - 11564.329(1) cm"' = 582.483(2) cm"'. However, no 
two measured states of ThO^ are separated by this energy. 
Because these states had not been observed by Goncharov et 
al. im, it is possible that the ThO^ population distribution in 
the discharge is superthermal, resulting in the observation of 
high-lying electronic states. Einally, we observed a ThO^ char¬ 
acteristic band head near 9800 cm"', but it was sufficiently out 
of range of our cw Ti;Sapph that we did not scan enough of the 
band to attempt a fit. 

4. Conclusion 

We have determined that X ^Ai is the electronic ground state 
of ThE^, and a is the first excited state, located 314 cm"' 
above the ground state. This was assigned by observing two 
Q = 0^ states that couple to both X and a, at T^. = 10,472 
and 14,296 cm"'. In an eEDM experiment it may be beneficial 
to first selectively ionize neutral ThE into the ThE^ state, 
then transfer population to the ground state to obtain a pure 
population in a single Zeeman, Stark, and hyperfine sublevel 
112 . Either Q = 0^ state can be used to transfer population 
between the two states. Transitions to the Q = 0^[14.29] state 
are nearly a factor of 10 stronger in relative intensity than to the 
Q = 0^[ 10.47] state; however, the weaker transition is notably 
on par with the transition strength to the ^no+ state observed in 
HfE^, which is already strong enough for coherent population 
transfer l2T|. 

Eor the eEDM experiment, the A-doubling in ^Ai is im¬ 
portant for setting the voltage required to fully polarize the 
molecule. The A-doubling in the ThE^ ^Ai state is about 7 
times larger than in the HfE^ ^Ai state li2TII . The larger value 
can be explained by the higher density of states at low energy. 
The ^Ai A-doubling can arise via L-uncoupling with a nearby 
^Ho state, which can in turn couple to a E-state via a spin-orbit 
interaction ioiiiii. Barker et al. observed a ^Hq state at E = 
3,395 cm"' IflTlI . three times closer in energy to ^Ai than it is 
in HfE^. This would in turn lower the energy of ^Hi, another 
potential interaction pathway, relative to HfE^. Additionally, it 
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is not unreasonable to expect a lower-lying S state to account 
for the remaining difference. 

In addition, we have measured and assigned rotational and 
vibrational constants of the ^A 2 and states, and measured 
rotational constants, vibrational constants, and electronic ener¬ 
gies for several intermediate electronic states. These values will 
provide a benchmark for ab initio calculations. Derived con¬ 
stants for ThF^ electronic states are given in Table A com¬ 
parison to other ThF^ experimental and theory results is given 
in Table Other than ambiguity in the ground state assign¬ 
ment, excellent agreement is seen between this work. Barker 
et al. C3, and theory calculations from Denis et al. OS and 
Skripnikov et al. M- 

Future directions for a ThF^ eEDM experiment will in¬ 
clude ionization spectroscopy of neutral ThF to determine ef¬ 
ficient ion creation pathways. The ThF dissociation energy 
is greater than the ionization energy iniiisa, which may al¬ 
low for greater flexibility in highly-efhcient ion creation. One 
prospect is to prepare the neutral molecule in a core nonpen¬ 
etrating high-1 Rydberg state, where it can vibrationally auto- 
ionize without perturbing the ion core ll42l l43l l44l . Addition¬ 
ally, resonance-enhanced multiphoton dissociation (REMPD) 
spectroscopy will be performed on ThE^. Laser-induced flu¬ 
orescence (LIE) spectroscopy had previously been used as the 
primary detection method on trapped HfF^ Il45l . but REMPD 
yielded a roughly 100-times improvement in count rate Bbl . 
We expect to use REMPD as our primary detection scheme in 
trapped ThP+. 

In conclusion, we have demonstrated that comb-vms is a 
powerful tool for high-resolution extensive survey spectroscopy 
of molecular ions, in this case ThE^. Our wavelength coverage 
spans the 696 - 900 nm region, covering 3,260 cm ' of con¬ 
tinuous spectra. Additional spectral regions in the visible and 
near-IR could be accessed with additional broadening and a dif¬ 
ferent VIPA etalon. 
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Table 2: Derived constants for observed states in ThF^ in cm ^ To is the origin of the v = 0 level relative to the v = 0 level 
of X ^A]. To errors are directly measured in transitions from ^A] and otherwise, error bars are adapted from El Vibrational 
constants are calculated assuming a Morse potential unless otherwise noted. Quoted errors are 95%. 


State 

To 

T 

A e 


a>e 

(^eXe 

a. [10-3] 


0 

0 

0.24311(7) 

656.96(1)+ 

1.920(3)+ 

1.00(4) 


314.282(7) 

314.0(2) 

0.24601(8) 

657.90(32) 

2.26(16) 

1 .12(6) 

^A2 

1052.5(1.0)1' 

1052.5(1.0)* 

0.24342(2) 

657.38(22) 

2.13(10) 

1.06(4) 

3A3 

3150(30)1' 

3149(30)* 

0.24388(7) 

659.31(29) 

2.00(15) 

1 .01(6) 

Q. = 0+[ 10.47] 

10471.889(7) 

10487.1(1) 

0.23583(5) 

626.67(56) 

1.88(17) 

0.97(7) 

ni[12.99] 

12994.3(1.0)* 

13032.6(1.0)* 

0.23049(2) 

580.33(10)° 

1.54(5) 

0.87(2) 

Hi [14.22] 

14223.6(1.0)* 

- 

0.23101(4)* 

- 

- 

- 

Q. = 0+[14.29] 

14295.967(8) 

- 

0.23187(5)* 

- 

- 

- 

Q.^0- 

14589.09(2) 

14620.81(1) 

0.22947(3) 

593.467(8)+ 

1.822(3)+ 

0.95(4) 

^n2 

16682(30)* 

16719(30)* 

0.22858(7) 

582.04(75) 

2.44(18) 

1 .22(6) 


Calculated without assuming a Morse potential. 

"Bq. 

^ Values and error bars from (m 

^ Calculated using ^A 2 and ^A 3 energies and errors from nn 
Morse potential assumption breaks down for v' = 2 and 3. 


Table 3: Comparison of low-lying states. Measured Tg, Bg, and coe values are com¬ 
pared from this work to measured values from Barker et al. El- and to theory calcula¬ 
tions from Denis et al. (m and Skripnikov et al. (m Quoted uncertainties ai‘e 95%. 


State 

Constant 

This Work 

El 

El 

El 

^Ai 

T, 

0 

315(1) 

0 

0 



0.24311(7) 

0.246(10) 

0.2439 

0.2438 


OJe 

656.96(1) 

658.3(2.0) 

667.3 

658.4 

'2+ 

Te 

314.0(2) 

0 

318.99 

448* 


B, 

0.24601(8) 

0.245(10) 

0.2464 

- 


OJe 

657.9(3) 

656.8(2.0) 

670.8 

- 

^A2 

Te 

— 

1052.5(1.0) 

1038.94 

1186* 


B, 

0.24342(2) 

0.243(10) 

0.2438* 

- 


OJe 

657.38(22) 

656.5(2.0) 

667.0* 

- 

3A3 

Te 

_ 

3150(30) 

3161.99 

3193* 


B, 

0.24388(7) 

- 

0.2443* 

- 


OJe 

659.33(26) 

665(30)* 

669.0* 

- 


* AGi/2. 

^ Stefan Knecht, unpublished results El 

^ Preliminary data for vertical transitions at = 3.75 a.u. El- The accuracy of 
the calculations has since been improved. 
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Table 4: Fitted constants for observed transitions in ThO^. Errors are statistical only and 
at the 95% confidence level. Bands were fit with half-integer J. No A-doubling terms were 
measured. Statistical errors on band origins are < 60 MHz, but an additional 150 MHz is added 
to include wavemeter error. No electronic states were able to be conclusively assigned. 



Vo 

B” 

B' 

D” 

[10-’] 

D' 

[lo-’i 

? (0, 0) 

12146.812(6) 

0.33984(2) 

0.32610(1) 

1.91(2) 

1.89(2) 

?(1,1) 

12091.106(7) 

0.33843(4) 

0.32469(4) 

1.78(7) 

1.77(6) 

?(1,0) 

13006.858(7) 

0.33973(4) 

0.32465(4) 

1 .6(1) 

1 .6(1) 

? (0, 0) 

11564.329(6) 

0.33900(2) 

0.32612(2) 

1 .86(2) 

1.87(2) 

?(1,1) 

11510.656(7) 

0.33767(3) 

0.32480(3) 

1.86(5) 

1.87(5) 
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